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Abstract 

In this paper, we focus on a special class of ideal projectors. With the aid 
of algebraic geometry, we prove that for this special class of ideal projectors, 
there exist "good" error formulas as defined by C. de Boor. Furthermore, we 
completely analyze the properties of the interpolation conditions matched 
by this special class of ideal projectors, and show that the ranges of this 
special class of ideal projectors are the minimal degree interpolation spaces 
with regard to their associated interpolation conditions. 
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1. Introduction 

The problem of polynomial interpolation is to construct a function p 
belonging to a finite-dimensional subspace of F[a:;] that agrees with another 
given function / on a set of interpolation conditions, where ¥[x\ := F[xi, . . . , x^] 
denotes the polynomial ring in d variables over the field F. If there exists 
a unique solution of the interpolation problem for every /, we say that the 
interpolation problem is poised. It's important to make the comment that F 
is a field of characteristic zero in this paper, for example F = Q, M, C. 

Error formulas for polynomial interpolation give explicit representations 
for the interpolation error. C. de Boor jsj derived a formula for the inter- 
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polation error. In terms of Newton fundamental polynomials^ T. Sauer and 



Y. Xu [ill, Il2|] presented Sauer- Xu error formulas for polynomial interpo- 
lation whose interpolation conditions have certain constraints. Afterward, 
C. de Boor 0] discussed the error formulas in tensor-product and Chung- 



Yao interpolation. In 1998, S. Waldron [18| investigated the error in linear 
interpolation at the vertices of a simplex. 

As an elegant form of multivariate approximate, ideal interpolation pro- 
vides a natural link between polynomial interpolation and algebraic geome- 
try. According to G. Birkhoff's definition [2|, a linear idempotent operator P 
on F[a?] is called an ideal projector if kerP is an ideal. In the theory of ideal 
interpolation, we are interested in finite-rank ideal projectors. A finite-rank 
ideal projector refers to the ideal projector whose range is a finite-dimensional 
subspace of ¥[x\. As mentioned by Carl de Boor, one reason for choosing 
ideal interpolation in the first place is the resulting possibility of writing the 
error formulas as in the following definition. 

Definition 1. (J^J) Let P be an ideal projector and {/ii, . . . , be an 
ideal basis for kerP. We say that the basis {hi, . . . , hm} admits a "good" 
error formula if there exist homogeneous polynomials Hj and linear operators 
Cj : F[a;] F[a;], j = 1, . . . , m such that for all / G ¥[x\, 



Hj{D)hk = Sj^k for all j,k = 1 

and 

m 

f-Pf = Y,C,{H,{D)f)h 



, m 



'31 



where Hj{D) will be defined in Section 2. 

We say that P has a "good" error formula if there exists an ideal basis 
{/ii, . . . , /im} for kerP that admits a "good" error formula. 

It is no surprise that every ideal projector in univariate polynomial ring 



has a "good" error formula [13|, [IJ] . When we turn to multivariate interpo- 
lation, things change greatly. C. de Boor [3] proved the existence of "good" 
error formulas for tensor-product and Chung- Yao interpolation. However, B. 



Shekhtman [15| showed that for a specific form of ideal interpolation by lin- 
ear polynomials in two variables, such a "good" error formula doesn't exist. 
Hence, the study of the type of ideal projectors with "good" error formulas 
is a rather complicated topic. 
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In this paper, we deal with a special class of ideal projectors, and prove 
the existence of "good" error formulas for this class of ideal projectors. It 
should be noted that the construction of the linear operators Cj,l < j < m, 
as in Definition [1] will be the subject of our future work. Moreover, we discuss 
the properties of the interpolation conditions matched by this special class 
of ideal projectors. The main results of this paper will be put in Section 3. 
The next section. Section 2, is devoted as a preparation for this paper. 

2. Preliminaries 

In this section, we will introduce some notation and recall some basic 
facts about ideal interpolation and algebraic geometry. For more details, we 
refer the reader to [s], Q, 0, El- 

Throughout the paper, we use N to stand for the set of nonnegative 
integers, and use boldface letters to express tuples and denote their entries by 
the same letter with subscripts, for example, ck = (ai, . . . , ad)- For arbitrary 
a. G N"', we define cxl = ail . . . a^!. 

A monomial x"" is a power product of the form x"^ . . . x'^'^ with ex. G N'^. 
We denote by T{x) := T(xi, . . . ,Xd) the set of all monomials in ¥[x]. For 
a polynomial f{x) = Cq^c" ^ ¥[x] with 7^ G F, we write the 

associated differential operator for / in the form 

f{D) = J2 

where 

aai-\ had 

= — 

dx'^' . . . dxf 

Henceforward, we use < to denote the usual product order on N*^. For 
a, (3 E N*^, a < (3 ii and only if < z = 1, . . . , d. In particular, a < /3 if 
and only if a < /3 and ex. ^ f3. A finite subset ^ C N'^ is lower if for every 
aeA, 0<l3<a implies /3 e A. 

A finite monomial set O C T{x) is called an order ideal if it is closed 
under monomial division, namely t & O and t'\t imply t' G O. For an order 
ideal O C T(a;), the corner set of O, denoted by C[0], is the set 

C[0] = {te T{x) : t ^ 0,Xi\t t/xi eO,l <i <d}. 
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Fix a monomial order -< on T{x), for all 7^ / G F[a3], we may write 

/= 0^(1)03' +C^(2)X^ H \-C^(r)X^ 

where ^ c^w G F,7W e i = 1, . . . , r, and a;T*" ^ £C^'" y ■ ■ ■ y x'^^^\ 
We shall call LT^(/) := c^(i)X'^''^^ the leading term and LM^(/) := a?'''*^' the 
leading monomial of /. 

Given an ideal I and a monomial order -<, there exists a unique reduced 
Grobner basis for X w.r.t. -<. Suppose that = {gi, . . . , gm}, then the 
set 

:= {a;" E T{x) : LT^{gj) f a;", for all 1 < j < m} 

is called the Grobner escalier of X w.r.t. -<. From the theory of Grobner 
bases, we know that M^iX) is an order ideal, and 

C[Ar^(X)] = {LT^((?i),...,LT^((7„)}. 

If P is a finite-rank ideal projector on F[a;], then there are two important 
subsets of F[a;] associated with P. The range of P is defined as 

V := ranP = {p E F[x] : p = Pf for some / G F[a3]}, 

which is a finite-dimensional subspace of F[a3], and the kernel space of P 

kerP = {g e ¥[x] : Pg = 0}, 

which forms a zero- dimensional ideal in F[a;]. Furthermore, as an infinite- 
dimensional F- vector space, F[a:;] has a corresponding dual space (F[a;])'. An 
ideal projector P on F[a;] also has a dual projector P* on (F[a;])', and the 
range of P* is 

A := ranP* = (kerP)^ = {A G (F[x])' : kerP C kerA}. 

Indeed, A is the set of interpolation conditions matched by P. It's easy to 
see that dim A = dim V and 

kerA := {/ G ¥[x] : A(/) = for all A G A} = kerP 

which satisfies 

kerA nV = {0}. 

The following theorems summarize some of the simple properties of ideal 
projectors. 
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Theorem 1. (f^J) A linear operator P -.Fix] — )> ¥[x] is an ideal projector if 
and only if the equality 

Pifg) = PifPg) 

holds for all f,gE ¥[x]. 



Theorem 2. ( ll d . \121) A linear operator P : ¥[x] — t- ¥[x] is an ideal pro- 



jector if and only if the operator P' := I — P satisfies 

P'ifg) = fP\9) + P'UPg) 

for all f,ge¥[x\. 
3. Main results 

In this section, we will describe a special class of ideal projectors with 
"good" error formulas in terms of ideal bases and interpolation conditions 
respectively. 

3.1. Representation in terms of ideal bases 



Following T. Sauer [10[, we refer to a reduced Grobner basis G for an ideal 
X as a universal Grobner basis if G is a unique reduced Grobner basis for 
X, independent of the monomial order. Now, we begin with an easy lemma 
about the universal Grobner bases. 

Lemma 3. // the ideal keiP has a reduced Grobner basis G = {gi, . . . ,gm} 
w.r.t. some monomial order, and the polynomials of G have the form 

Qj = ^"'^^ — ^ Cj,/3CC^, for 1 < j <m and Cj^f} G F, (1) 

0</3<aO) 

then G is a universal reduced Grobner basis for keiP w.r.t. any monomial 
order, and the monomial set 

= {x'^ e T{x) : x"''' t £c^, for all 1 < j < m} (2) 

is the unique Grobner escalier o/kerP w.r.t. any monomial order. 
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Proof. For an arbitrary j with 1 < j < m and an arbitrary /3 G N'^ with 

< /3 < ct^^\ we have that \ x°'^^\ Suppose that -< is an arbitrary 
monomial order, then x^^ \ x""^^^ together with f3 7^ a^-'-* imphes x^ -< a:;"'^ . 
Consequently, for an arbitrary monomial order -<, l!T^{gj) = a;"*^^ with 

1 < j < m, and ^-polynomial of Qi and gj with 1 < i < j < m is the 
combination 

S{9r,9j) = 9^ 9J, 

where LCM{x°'^'\ x"^^^) is the least common multiple of x"^'^ and 

Since G is a reduced Grobner basis for kerP w.r.t. some monomial or- 
der, it follows that S{gi,gj) reduces to zero module G w.r.t. this monomial 
order. Indeed, for arbitrary monomial order -<, LT^{gj) = x"-^^\ it implies 
that S{gi,gj) reduces to zero module G w.r.t. any monomial order. There- 
fore, we can say that G is a universal reduced Grobner basis for kerP w.r.t. 
any monomial order. Furthermore, it follows that O is the unique Grobner 
escalier of kerP w.r.t. any monomial order. □ 

Proposition 4. Let {gi, . . . ,gm} be a reduced Grobner basis for kerP sat- 
isfying condition (Cp, and O be a monomial set as in Then for every 
monomial x'^ G T{x), there exist polynomials A~yj, 1 < j < m such that 

m 

P'{x'^) = x'^-P{x'^) = Y,A^,m (3) 



and 

In other words, 



A^j = if aj^-'^fx-^. (4) 



P'(x'^) = - P(x^) = Af,9r (5) 

Ql(j)<7 

Proof. For every 7 G N'^, define an ideal 

= {9j ■■ < 7) 

To prove this proposition, it suffices to show that x'' — P{x'^) G J-y for every 
x'^ G T[x). Assume not and let x'^ be a monomial of least total degree such 
that x"^ — P{x'^) ^ J-y. Since for every x^ G O, 

= £c^ - P(a;^) G J^, 
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we know that x"^ ^ O. Therefore, we can find some 1 < j < m such that 
QjO) < ^ Let (5 = 7 — cx^'^^ > 0. By Lemma El we have 



x^gj = x^ - ^ Cj^fiX' 



Consequently, 



x'g, = P\x'g,) = P'{x^) - Yl e J^- 

But for every f3 such that f3 E O, f3 < a^^^ we have /3 + 5 < 7- Recall 
that a;''' is a monomial of least total degree such that P'{x'^) ^ J~^. Hence, 
P\x^^^) G Jfi+s C Xy. Since X,, is an ideal, then P\x'') G Xy. This is a 
contradiction to our hypothesis. □ 

We need a standard key lemma for factorization of homomorphisms. 



Lemma 5. n i j . uq/) Lei A : X ^ Y and B : X ^ Z be two linear operators 
between linear spaces X,Y and Z. Then there exists linear operator C such 
that 

A = CB 

if and only if 

keiB C kerA. 

The fact that an ideal projector P has a "good" error formula depends 
on not only the ideal basis for kerP, but also the choice of ranP. Next is the 
main theorem of this paper, which states that the ideal projectors satisfying 
the conditions of Theorem O have "good" error formulas. 

Theorem 6. Suppose that an ideal kerP has a universal reduced Grobner 
basis G satisfying the conditions of Lemma 0, and ranP is 

V = span^{x'^ G T{x) : £c"''^ f x^, for all 1 <j < m}. (6) 

Then G is the ideal basis for kerP that admits a "good" error formula. 

Proof. Define operators Aj on T{x) by letting 

( A ■ if rv^-?) < -y- 
A^{x'^) = \ (7) 
(J, otherwise. 
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where Aj^^ are defined in and extend Aj^-^ by linearity on ¥[x]. Then by 
dS]) and linearity, we have 

f-Pf = P'f = Aj{f)gj. 

By dZD, 

keiAj D span{a;^ : a;"'" f £c^} = ker(-j-yyD°^'"). 

Hence, by Lemma [5l there exist operators Cj such that Aj = Cj o Hj{D) 
where Hj{D) := ^^D"'^'. It is trivial to check that Hj{D){gk) = 5j^k- □ 

In the following, we will present some examples to illustrate the conclusion 
of Theorem [61 

Example 1. Let P be a Lagrange projector onto spauipjl, xi, X2, Xg} with 
the interpolation point set {(1, 0), (1, 1), (1, 2), (2, 0)} C F^. Then by Theo- 
rem [6l 

{(Xi - - 2),X2(X2 - 1)(X2 - 2),X2{Xi - 1)} 

is the ideal basis for kerP that admits a "good" error formula. □ 

Example 2. Let P be an ideal projector onto spaupjl, xi, X2, xf, X1X2, x\, x\} 
given by 

Pxlx2 = 

Pxl = X2 

Pxix\ = X\X2 
TDrfA _ 9^3 _ 2 

Then the ideal basis 

r2 23 32 24 4~1 

admits a "good" error formula. □ 

Example 3. Let P be a Lagrange projector onto the spauj^jl, xi, X2, X3} 
with the interpolation point set {(0, 0, 0), (0, 1, 0), (0, 0, 1), (1, 0, 1)} C 
Then 

{xiX2, X2X3, xj - Xi, x\ - X2, X3 - X3, X1X3 - Xi} 

is the ideal basis for kerP that admits a "good" error formula. □ 
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We select test functions 

/l(xi,X2) = (l-Xi)2 + (l-X2)' + l, 

/2(Xi,X2) = xl+xl, 
h{xi,X2, X3) = (1 - Xif + (1 - X2f + (1 - 0:3)^ + 1, 
/4(xi, X2, X3) =xl+xl + xl 

to illustrate the "good" error formulas about the ideal projectors in the above 
examples. 

For Example 1, we have 

/i-P/i = (xi-l)(xi-2), 

/2 - Pf2 = (xi + 3)(a:i - l)(xi - 2) + X2(x2 - l)(x2 - 2). 
For Example 2, we get 

/i - Pfi = 0, 

/2 - Pf2 =xl- X2. 

For Example 3, 

/a - Pf'i = xl- xi+ xl- X2 + xl- X3, 

/4 - P/4 = (xi + l){xl - Xl) + {X2 + - 0:2) + (a;3 + l){xl - x^). 

3.2. Representation in terms of interpolation conditions 

Next, we will describe the properties of the interpolation conditions matched 
by the ideal projectors satisfying the conditions of Theorem [6l 

Proposition 7. Let P he an ideal projector, and A = {Ai, . . . , A„} the set 

of interpolation conditions matched by P. Let -<iex{i), 1 < i < d, be the 
lexicographic order 

Xi >-■■■>- Xd y Xl y ■■■ >- Xi^i. 

Then kerP has a universal reduced Grobner basis G satisfying the conditions 
of Lemma if and only if 

A/'-^,e.(i,(kerA), Ar,,,(,,(kerA), . . . , Ar,,^^,,,(kerA) 

are identical. 
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Proof. One direction of the proof is obvious due to the fact that kerP has 
a universal reduced Grobner basis G satisfying the conditions of Lemma |31 

To prove the converse, assume that G^^^^^.^, 1 < z < d, is the reduced 
Grobner basis for kerA w.r.t. -<iex{i)- Indeed, if 

A/'-...(.)(kerA) = AT,,^^,,, (kerA) = ■ ■ ■ = AT,,^^^,, (kerA) = O, 

then it is easy to prove 

Suppose that G = {gi, . . . ,gm} and C[0] = {x"*^', . . . , a;"'"'}. Then rear- 
ranging the elements of C[0] appropriately, we have 



^T^^^xaM) = ^T^^.xM = ■■■ = LT^,,,< J^7,) = VI < j < 



m. 



Since for arbitrary fixed 1 < z < d, G is the reduced Grobner basis for kerA 
w.r.t. -<iex{i)i it follows that the polynomials of G have the form 

Qj = a;"'^* — Cj^fjX^, for all 1 < j < m. 

.0) 



<lex(i)^ 



Furthermore, by the property of lexicographic order, for arbitrary 1 <i < d, 
x^~<iex{i)X°'^^^ implies that < (3 < cxS^\ Hence, we can deduce that the 
polynomials of G have the form: 

Qj = a;"*^'' — Cj^fiX^ , for all 1 < j < m. 

0</3<q;(j) 

This completes the proof. □ 

Moreover, Proposition [7] coupled with Theorem |6] immediately implies the 
following useful corollary. 

Corollary 8. Let P he an ideal projector, and A = {Ai, . . . , A„} the set of 

interpolation conditions matched by P. If 

(kerA) = Af,^^^^^^ (kerA) = ■ ■ ■ = Af,^^^^^^ (kerA) = O, (8) 

where A/'^,^^(;j(kerA), 1 < i < d are as above, then the ideal projector P onto 

V = spanp{x^ ■.x'^eO} (9) 
has a "good" error formula. 
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Remark 1. Let . . . , ^(^^^ E be distinct points and A^'^\ A'^f'^ C N'^ 
lower sets. Suppose that the set of interpolation conditions has the form 

A = {6^w o : a e for all 1 < A; < //}, 

where S^^k) denotes the evaluation functional at the site ^('^^ , then J\f^^^_^^.^ (kerA) 
can be directly computed by the fast algorithms given in p, l8| without com- 
puting the Grobner basis for kerA. 

Example 4. Suppose that the set of interpolation conditions matched by P 
is as follows: 

A = {5(0,0), Vo) ° D^'''\ ^0,0) o D^'^'\ 6^0,1), ° ^^'""^ '^(i,o), '^(i.o) o D^'''^}. 
Since 

(kerA) = Af^^^^^^.^(keTA) = {1, X2, xj, xi, xiX2, xj, x^}, 
then P onto spanjpjl a "good" error formula. □ 

As mentioned by Carl de Boor in j^, the existence of a "good" error 
formula for an ideal projector restricts the range of ideal projector to be of 
least degree. The following theorem is a particular case of this fact. Here, 
we also provide a simple proof, for completeness. 

We denote by 11^ the subspace of polynomials in ¥[x] of total degree 
at most r. Suppose A = {Ai,...,A„} C and / G ¥[x], we write 

A(/) = (Ai/, A2/, . . . , Xnff. For a finite set F = {/i, ■ ■ ■ , /,} C ¥[x], A(F) 
signifies the n x k matrix whose columns are A(/j), 1 < i < k. 

Theorem 9. Let A = {Ai,...,A„} be the set of interpolation conditions, 
If A satisfies condition then V defined in ^ is the minimal degree 
interpolation space w.r.t. A. 

Proof. Suppose that the maximal total degree of the monomials in O is r, 
then V C n^. It's obvious that the interpolation problem of finding p E V 
such that 

XiP = Xif, for all 1 < z < n 

is poised. According to T. Sauer's definition (cf. j^), we need to prove two 
properties on this special class of projectors. Firstly, the operator P onto V 
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is degree-reducing, namely for each / G 11^ with < k < r, the interpolating 
polynomial Pf also belongs to 11^. Secondly, the subspace C 11^ is of 
minimal degree, namely there is no subspace V C Ilr-i such that the above 
interpolation problem is poised. 

Since each / G Ilfc can be written in the form 



then 

Since V is the range of P, we have that for any G O, 

P(x^) = x^. (11) 

On the other hand, if x^ ^ O, then there must exist some a;" G C[0\ 
such that x"-\x^. From Corollary |8l it follows that for some c-,, G F with 
< 7 < a, 

— c-^x'' G kerP. 

0<7<q; 

Multiplying the above equation by a3^~", we get 



X 



,/3 



c^xi'+^-° G kerP 



/3-Q;<7+/3-a</3 



If x''^^ " ^ (9, we repeat the above processing. Finally, we can find some 
x^ ^ O with (3' < (3, and associated coefficients Cf^i G F such that 

a;^ - ^ Ci^ix'^' G kerP. 



/3'</3 



Since a;'^ G Ilfc, it follows that 



P(a;^) = ^ nfc. (12) 



/3'<,9 
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From fllUI) . fill I) . fll2p . we can conclude that for arbitrary / G 11^ with 1 < 
k<r,Pfe Hfc. 

To prove the minimal degree property, we set O' = 0[ [JO2, where 
0[ := {x^ e T{x) : x^ E Ur^i and x^ E C}, 

and 

O'^ := {x^ E T{x) : x^ E n^_i and x^ ^ O}. 

Then we need only to prove that the matrix A{0') has rank less than n. 

Recalling equality ( IT^ . we can easily see that for an arbitrary x^ E O2, 
A[x^) linearly depends on the columns of A{0[). Equivalently, A{0') has 
rank less than or equal than ^0[. Since at least one x^ E O belongs to 11^ 
and not to n^.i, we have that the matrix A[0') has rank less than n. To 
sum up, we can say that V = span^ja;^ : x^ E O} is the minimal degree 
interpolation space w.r.t. A. □ 

Acknowledgements 

The authors wish to thank the anonymous reviewer for valuable sugges- 
tions and comments that have improved the presentation of the paper. 

References 

[1] T. Becker, V. Weispfenning, Grobner Bases, Vol. 141 of Graduate Texts 
in Mathematics, Springer- Verlag, New York, 1993. 

[2] G. Birkhoff, The algebra of multivariate interpolation, in: C. V. Coff- 
man, G. J. Fix (Eds.), Constructive Approaches to Mathematical Mod- 
els, Academic Press, New York, 1979, pp. 345-363. 

[3] C. de Boor, On the error in multivariate polynomial interpolation, Appl. 
Numer. Math. 10 (1992) 297-305. 

[4] C. de Boor, The error in polynomial tensor-product, and Chung- Yao, 
interpolation, in: A. LeMehaute, C. Rabut, L. Schumaker (Eds.), Sur- 
face Fitting and Multiresolution Methods, Vanderbilt University Press, 
Nashville TN, 1997, pp. 35-50. 



13 



C. de Boor, Ideal interpolation, in: C. K. Chui, M. Neamtu, L. L. Schu- 
maker (Eds.), Approximation Theory XI: Gatlinburg 2004, Nashboro 
Press, Brentwood TN, 2005, pp. 59-91. 

L. Cerlienco, M. Mureddu, From algebraic sets to monomial linear bases 
by means of combinatorial algorithms. Disc. Math. 139 (1-3) (1995) 73- 
87. 

D. Cox, J. Little, D. O'Shea, Ideal, Varieties, and Algorithms, third ed., 
Undergraduate Texts in Mathematics, Springer, New York, 2007. 

B. Felszeghy, B. Rath, L. Ronyai, The lex game and some applications, 
J. Symb. Comput. 41 (6) (2006) 663-681. 

T. Sauer, Polynomial interpolation of minimal degree, Numer. Math. 
78 (1)(1997) 59-85. 

T. Sauer, Lagrange interpolation on subgrids of tensor product grids. 
Math. Comp. 73 (245) (2004) 181-190. 

T. Sauer, Y. Xu, On multivariate Lagrange interpolation. Math. Comp. 
64 (1995) 1147-1170. 

T. Sauer, Y. Xu, On multivariate Hermite interpolation. Adv. Comput. 
Math. 4(1995) 207-259. 

B. Shekhtman, On one question of Ed Saff, Elec. Trans. Numer. Anal. 
25 (2006) 439-445. 

B. Shekhtman, Ideal interpolation: Translations to and from Alge- 
braic Geometry, in: L. Robbiano, J. Abbott (Eds.), Approximate Com- 
mutative Algebra, Texts and Monographs in Symbolic Computation, 
Springer- Vienna, New York, 2009, pp. 163-192. 

B. Shekhtman, On non-existence of certain error formulas for ideal in- 
terpolation, J. Approx. Theory 162 (7) (2010) 1398-1406. 

B. Shekhtman, On the naive error formula for bivariate linear interpo- 
lation, in: Wavelets and Splines: Athens 2005, Mod. Methods Math., 
Nashboro Press, Brentwood, TN, 2006, pp. 416-427. 



14 



[17] B. Shekhtman, On error formulas for multivariate polynomial interpola- 
tion, in: M. Neamtu, L. Scliumaker (Eds.), Approximation Theory XII: 
San Antonio 2007, Nashboro Press, Brentwood TN, 2008, pp. 386-397. 

[18] S. Waldron, The error in linear interpolation at the vertices of a simplex, 
SIAM J. Numer. Anal. 35 (3) (1998) 1191-1200. 



15 



